On the temperature dependence of the superconducting gap in high-T c cuprates 
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It is proposed that (i) the temperature dependence of the superconducting gap A(T) in high-T c 
cuprates can be predicted just from the knowledge of A(0) and the critical temperature T c ; and, in 
particular, (ii) A(0)/T c > 4 implies that A(T C ) / 0, while A(0)/T c < 4 corresponds to A(T C ) = 0. A 
number of tunneling experiments appear to support the above proposition, and, furthermore, show 
reasonable quantitative agreement with a model based on the two-dimensional stripe hypothesis. 

PACS numbers: 74.50.+r, 74.72.-h, 74.20.Mn, 74.25. Jb 
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In conventional superconductors, the phenomenology 
of the superconducting (SC) energy gap closely follows 
that of the model of Bardeen, Cooper and Schrieffer 
(BCS)Q. Namely: ( i) Below the SC transition tempera- 
ture (T c ), there are no one particle excitations within the 
SC gap A(T), where T is the temperature; (ii) At the 
gap boundaries ±A(T), there are symmetric SC peaks; 
(hi) The ratio A(0)/T c is close to 1.76; (iv) A(T C ) = 0. 
In contrast, the tunneling studies of high-T c cuprates 
show that: (i) The density of states within the SC gap 
is not zero; (ii) The SC peaks are frequently asymmet- 
ric; (iii) The ratio A(0)/T c is significantly larger than 
1.76 and not universal; (iv) There are indications that, 
in some high-T c compounds, A(T C ) ^ 0. The above non- 
universality implies that any BCS-like one-parameter de- 
scription of the tunneling spectra cannot be adequate for 
high-T c cuprates. The complexity of these materials fur- 
ther suggests that their proper microscopic description 
should contain many parameters. In this work, however, 
I report a surprising finding that the phenomenology of 
the SC gap in high-T c cuprates may be describable by 
only two parameters. 

The present study was originally undertaken to test 
a simple model proposed in Ref. |2|. The calculations 
of the tunneling characteristics based on that model re- 
quire only two input parameters A(0) and T c . The rea- 
sonable success of this test, of course, lends credibility to 
the model, but also, irrespectively of the model, it sug- 
gests that, even though the ratio A(0)/T c may exhibit 
a non-universal doping dependence, and, moreover, may 
be different for different samples of the same material, 
the measured values of A(0) and T c are sufficient to pre- 
dict the entire evolution of A(T). Of particular interest 
is the following model-based rule: 

A(0)/T c > 4 & A(T C ) ^ / asymmetric SC peaks; 
A(0)/T c < 4 <S> A(T C ) = / symmetric SC peaks. 

(1) 

It was the surprisingly good agreement of the 
symmetry /asymmetry part of this rule with many 
experiments |2j, that prompted me to test the model pre- 
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FIG. 1: 2D configuration of diagonal stripes (a cartoon) 



dictions for the shape of the A(T)-curve. 

The model of Ref. Q is based on the hypothesis that 
a two-dimensional (2D) arrangement of stripes shown in 
Fig- d exists in the Cu02 planes of high-Tc cuprates. In 
such a superstructure, there exist hole states localized ei- 
ther inside the stripes (b-states) or in the antiferromag- 
netic (AF) domains between the stripes (a-states). Due 
to the small size of the AF domains and stripe elements 
(20-30A), the energy levels within each unit should be 
broadly spaced (estimated 2] as 40 meV). The model pro- 
ceeds by assuming the existence of only one low-energy 
a-state per AF domain and two low-energy b-states per 
stripe element. As argued in Ref. Q, a promising choice 
for the model Hamiltonian is 

)7<=1 

+ 9 Y.( b t+ b t- a > a J+ h - c -)' ( 2 ) 
where index i or j labels AF domains; implies, that 
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the jih AF domain is the nearest neighbor of the ith 
domain; a; is the annihilation operator of a hole inside 
the ith AF domain; 6^ CT is the annihilation operator of a 
hole inside the stripe element separating the ith and the 
jth AF domains; a is the spin index, which can have two 
values "+" or "— "; e a and Eb are position-independent 
on-site energies for a- and b-states, respectively, counted 
from the chemical potential; and, finally, g is the coupling 
constant. The spins of a-states are tracked by index rji, 
which alternates between 1 or —1 thus tracing the sign of 
the AF order parameter in the AF domains. The sum su- 
perscript "r/i = 1" indicates that the summation extends 
only over the supercells having r\i = 1. 

The total energy of such a system as a function of the 
chemical potential fi has at least two minima^: H = £ a 
and /i = £t,, which implies two most promising regimes 
- Case I: Eb = 0; and Case II: e a = 0. The mean-field 
solution then leads to the density of states described by 
the following system of equations^: in Case I, 
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e j4 (k) = \\ei 



g 2 (2n B - 1) 



1 



g 2 (2n B -l) 2 \V(k)\ 2 , 



(3) 



64tt 2 
and. in Case II, 



7T p 7T 



7T J — 7T 



dk v I dk, 



(l-2ru(k)|V(k)| 5 



6A(k) 



Mk) = - 



g 2 \V(k)\ C a (l-2n B ) 
8e B 



£B = ^ + 16 g2 °« 



(4) 

(5) 
(6) 



where subscripts A and B indicate the Bogoliubov quasi- 
particles associated with a- and b-states, respectively; 
£ a/b and jia/b are the quasiparticle energies and the 
occupation numbers, respectively. A-quasiparticles are 
characterized by well-defined wave vectors k. All B- 
quasiparticles have the same e B - The expressions for 
n J 4(k) and n B , and for the auxiliary quantities C a and 
V(k) are: 



ru(k) 



1 



exp 



n B 



exp(^) + 1' 



V(k) = 2 



(7) 



(8) 



(9) 



Ca = i/ dkx j dk v( 2n ^) - !)l^( k )l- (10) 



The tunneling spectra of A-states have Van Hove singu- 
larities, which, in Case I, are located at 



£j40 = ± y/el+g 2 (2n B -l) 



and, in Case II. at 



£.40 = ± 



g 2 C a (l-2n B ) 
<±e B 



(11) 



(12) 



Three representative zero-temperature tunneling spectra 
of both A- and B-states can be found in Ref . |2| . 

In both Cases I and II, the mean-field T c can be ob- 
tained ,2j numerically from the following equation: 
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exp 
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(13) 



Given the above picture, the basic features of the 
tunneling phenomenology can be interpreted as follows: 
(1) The SC peaks at ±A, are identified with the Van Hove 
singularities at ±eao in the density of A-states. (2) The 
density of B-states is assumed to be more difficult to de- 
tect due to, perhaps, smaller tunneling elements. 

The knowledge of A(0) and T c is thus sufficient to 
obtain the model parameters \e a \ and g in Case I 
[Eas. (|llll3|) with n B = 1], and |e;,| and g in Case II 
rEas. fciWIlOllfllfl with n A (k) = l;n B = 0]. After 
that, A(T) can be calculated numerically using Eqs.(|3- 
1121 without additional adjustable parameters. Impor- 
tantly, Cases I and II do not overlap in terms of the 
ratio A(0)/T c : in Case I, A(0)/T c > 4, while in Case II, 
A(0)/T c < 4. The entire family of the A(T)-curves is 
presented in Fig. [21 

In both Cases I and II, \e a — £&| is a measure of the 
pseudogap (see Sec.V of Ref. 0). This pseudogap is a 
real feature in the normal density of states (as, e.g., in 
the proposal of Ref.|J), i.e., it is not associated with the 
SC fluctuations above T c . This feature is a generic con- 
sequence of the nanoscale phase separation — two differ- 
ent kinds of real space regions (stripes and AF domains) 
imply two different energies associated with them. The 
model pseudogap equals zero only in the critical case, 
e a = Eb, which, as argued in Ref.|2|, may correspond to 
the doping concentration 0.19. The model has two differ- 
ent SC energy gaps: eao f° r A-states and e b for B-states. 
In Case I, e b turns to zero at T = T c , while eao becomes 
equal to the pseudogap \e a — e&|. In Case II, the situa- 
tion at T = T c is opposite: eao — 0, while e b = \e a — £&|- 
Thus the difference between the two model regimes is 
not qualitative but only quantitative. However, since, by 
assumption, the observable SC peaks represent only A- 
states, the phenomenology of the SC gap appears to be 
qualitatively different in the two regimes. 

In Fig. |3 I compare the theoretical predictions with 
the results of the break junction (BJ)0, 0, Q and the 
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FIG. 2: (Color online) Family of theoretical curves for the 
temperature dependence of the superconducting gap. Thick 
line corresponds to the critical ratio A(0)/T c = 4. Solid lines 
above the thick line describe Case I, and below the thick line 
Case II. The dashed line represents the BCS result. 



interlayer tunneling (ILT)@, E H El El E3 experi- 
ments. The ILT technique gives better energy resolu- 
tion and comes closer to measuring the true intrinsic 
density of states in the bulk of cuprates, but, at the 
same time, it generates significant overheating of the 
sample. In order to avoid this problem, the authors of 
Refs.0, E3 

have used a short-pulse method, which, 
however, comes at a cost of somewhat lower energy reso- 
lution [Fi gs-ETc-g) a nd (q-s)]. Alternatively, the authors 
ofRefs.||lirJjii|l3] [Figs.|lh-n) and (t)] have perfected 
the "conventional" ILT technique by manufacturing suf- 
ficiently small mesas and thus reducing the overheating. 
In the latter case, the experimental data points reflect the 
ambient temperature of experiment. The question, how- 
ever, remains concer ning the true value of the tempera- 
ture of the mesasfLU Il5|. Theoretical and experimental 
estimates made in Refs.[l3, llfl LLZ|> have indicated pos- 
sible significant overheating, but there is still a chance 
that, at least in the underdoped Bi-2212, the tempera- 
ture was not too much distorted at the voltage biases 
probing the SC peaks. 

Most of the theoretical curves shown in Fig. E] agree 
with the experimental data within the limits of [some- 
times, large] uncertainty associated with the broadening 
of the SC peaks (see the original references). The data 
points in Fig. E] are obtained from the peak-to-peak sep- 
aration in the experimental tunneling spectra. Usually, 
the peak broadening shifts the peak maximum towards 
higher energies. If this effect is undone by extracting 
the "bare" values of A, then the experimental points in 
Figs. E{a-g,o, q-s) would appear significantly closer to the 
theoretical curves. This is a consequence of the fact that 



the SC peaks are broader at higher temperatures. The 
magnitude of the above correction can be judged on the 
basis of the estimates of bare A made in Ref. and pre- 
sented in Fig.|3Ia,b). [Note: the widths of the SC peaks 
corresponding to Figs.E{a,b) are comparable or smaller 
than those corresponding to Figs.E{d-g,o-s).] 

The agreement between the theory and the experiment 
is typically better, whenever the SC peaks are sharper. 
This is, particularly, true for the "pulsed" ILT experi- 
ments (Fig. E{c-g) and (q-s)): up to T/T c ~ 0.7 the data 
points from the overdoped samples [Figs. EJc,g,s)] cor- 
respond to sharper SC peaks and agree better with the 
theory. Figure Efn) constitutes an exception from the 
above trend, but, in this case, a stronger degree of over- 
heating is suspected. The best agreement between the 
theory and the experiment can be observed in Figs. Ef h- 
m) , in which case the data points were generated by the 
"conventional" ILT technique (expected to give the best 
energy resolution). If these data also stand the test of 
time with respect to the true temperature of the sam- 
ples, then such an agreement would amount to a striking 
success of the present model. 

The experimental data in Figs. 01 are also fairly con- 
sistent with the first line of the rule 0J. As far as the 
second line of that rule is concerned, then Figs. E{i) and 
(n) can be cited as supporting the rule. The data points 
in Figs. Efc,g,o,s) while corresponding to A(0)/T c < 4 
show trend towards A(T C ) ^ 0. However this trend is due 
to the highest temperature points, which originate from 
poorly characterized SC peaks. The conclusive verifica- 
tion of rule Q should thus await further experiments. 
It should also be noted in this context, that the scan- 
ning tunneling spectroscopy studies of Bi-2212[18j and 
Bi-22010, which indicated that A(T C ) ^ 0, were all 
performed on the samples exhibiting A(0)/T c > 4. 

In a future study, I plan to test the model predictions 
for the specific heat experiments, which should be equally 
sensitive to the densities of both A- and B-states. 

In conclusion, I express cautious optimism that 
(i) there exists a two-parametric description of the A (In- 
curves in high-T c cuprates; (ii) there exists a critical ratio 
A(0)/T c , which signifies a transition between two differ- 
ent kinds of SC states, and (iii) the above critical ratio 
is approximately equal to 4. The quantitative agreement 
between the experimental data reproduced in Fig.EIh-m) 
and the model calculations also appears quite promising. 

The author is grateful to V. M. Krasnov, A. Yurgens, 
N. Miyakawa, M. Suzuki, and V. N. Zavaritsky for valu- 
able discussions. 
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FIG. 3: (Color online) Temperature evolution of the SC gap. Solid lines - theory; circles - experiments; BJ - break junction; 
ILT - interlayer tunneling. Horizontal marks in each frame indicate A/T c — 4. Whenever the original tunneling spectra were 
reported (a-h, k, n-t), filled circles indicate reasonably well-pronounced SC peaks, while open circles indicate very broad and 
small SC peaks. Filled circles in frames (i,j,l,m) imply no additional information about the SC peaks. Squares in frames (a,b) 
represent the "bare" values of A estimated in Ref.Q- 
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